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1.0  INTRODUCTION 


The  classical  transportation  problem  is  well  known  for  its  widespread 
applicability  and  for  the  facility  with  which  it  can  be  solved.  Many  pro¬ 
cedures  have  been  developed  for  reformulating  disparate  linear  progranraing 
problems  as  transportation  problans  to  take  advantage  of  the  computational 
efficiency  inherent  in  the  specialized  transportation  algorithm  [2,3,8,10]. 

We  generalize  the  concepts  of  an  earlier  paper  [6]  to  shew  how  constrained 
transportation  problems  (i.e. ,  transportation  problems  with  additional 
linear  constraints)  can  be  transformed  into  pure  transportation  problems. 

Our  procedure  determines  if  an  arbitrary  extra  linear  constraint  can  be  ' 
transformed  into  an  equivalent  bounded  partial  sum  of  variables  involving  a 
single  node  constraint.  If  this  is  possible  the  procedure  gives  the  linear 
transformation  that  yields  the  equivalent  constraint.  'Ihis  extends  the  work  of 
Wagner  [9],  Manne  [4,p.  382],  and  Charnes  [1]  who  have  shown  how  transportation 
problems  with  these  bounded  partial  sums  can  be  reformulated  into  pure  trans¬ 
portation  problems 

It  is  conjectured  that  our  transformation  requires  oonputational  effort  on 
the  same  order  as  that  required  to  find  an  initial  basic  solution  by  Vogel's 
Approximation  Method.  Since  ccmputational  results  indicate  that  specialized 
transportation  codes  can  solve  transportation  problems  at  least  150  times 
faster  than  general  purpose  linear  programming  codes  [5] ,  our  results  make  it 
possible  to  solve  constrained  transportation  problems  of  the  specified  class  with 
substantially  greater  efficiency  than  by  a  general  purpose  algorithm. 

2.0  PROBLEM  STATEMENT,  MATHEMATICAL  DEVELOPMENT,  AND  AN  EXAMPLE 


A  transportation  problan  with  one  additional  extra  constraint  (  a 
singularly  constrained  transportation  problem)  can  be  stated  mathematically 


-2- 


as  follows: 


Minimize 

2 

1t-M 

£  c1j  xij 

jcN  J  J 

subject  to 

I 

jeN 

»ij  ■  »i.  ><* 
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i  eM 

*ij  -  ^ 
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icM 

jcN  '  d 

Xjj  >  0  ,  1t.M, 
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..  ,m}, 

N  r[l,2,...n},  z 

kM 

jcN  J 

form  the  singularly  constrained  transportation  problem  appears  as  indicated  in 
Table  1.  The  coefficients  (p-jj)  of  the  extra  constraint  appear  below  the  familiar 
echelon-d1ago;ial  structure  of  ones  of  the  transportation  problem. 
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Table  1.  Oetached  Coefficient  Form  of  Singularly 
Constrained  Transportation  Problems 
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Our  goal  is  to  specify  a  computationally  simple  procedure  for  identifying  a 

linear  combination  of  the  ordinary  transportation  constraints,  if  one  exists,  which 
can  be  subtracted  fran  the  extra  constraint  to  produce  an  inequality  that  is  equi¬ 
valent  to  establishing  a  bound  for  a  partial  sun  of  variables  associated  with  a 
single  origin  or  destination.  Fran  well  known  properties  of  the  transportation 
matrix,  if  there  exists  a  linear  canbination  of  the  transportation  constraints 
that  has  the  desired  form,  then  a  linear  combination  can  be  found  in  which  any 
particular  transportation  constraint  receives  a  zero  weight.  Thus,  we  may  arbi¬ 
trarily  delete  a  node  constraint  when  seeking  such  a  linear  combination.  Having 
done  this  certain  variables  can  be  viewed  as  having  only  one  entry  in  the  coefficient 
matrix  for  purposes  of  finding  the  desired  linear  combination.  To  exploit  these 
facts  in  seeking  a  transformation  of  the  singularly  constrained  transportation  problem 
into  an  ordinary  transportation  problem,  we  arbitrarily  emit  the  first  origin  con¬ 
straints.  Our  principal  observations  are  then  contained  in  the  following  three  cases. 

Case  a.  Assume  that  the  extra  constraint  is  equivalent  to  a  partial  sun  of 
variables  in  a  siiigle  origin  constraint  q  other  than  origin  1  (i.e.,  the  equivalent 
extra  constraint  is  of  the  form  Z  x  .=f  for  Scn  as  illustrated  in  Table  1). 

Having  set  equal  to  0,  the  unique  values  for  all  the  destination  constraint  mul¬ 
tipliers  can  be  determined  using  the  equations  R}+Kj=p^j .  This  is  possible  by  the 
assumption  that  the  equivalent  constraint  invoices  variables  in  origin  q^l  (i.e. , 
the  coefficient  on  the  variables  ,  j=l,...,  n  in  the  equivalent  constraint  must 
equal  0).  By  the  same  reasoning  a  unique  origin  multiplier  can  be  found  for  every 
other  origin  except  origin  q.  In  origin  q  not  all  of  the  equations  Rq  +  Kj  = 
can  simultaneously  be  satisfied  by  a  single  value  for  There  will  be  two  values, 
one  that  satisfies  the  equation  for  a  subset  Sc  N  and  another  that  satisfies  that 
equation  for  the  destinations  in  the  subset  N-S.  By  setting  equal  to  the 
value  Rq  +  Kj  =  pqj  for  jf.N-S  a  linear  combination  of  the  standard  node  constraints 
has  teen  found  which  when  subtracted  ficm  the  original  extra  constraint  yields 
a  restriction  on  the  partial  sum  of  variables  in  origin  constraint  q  associated 


d-  i  R^a^  -  i  K_:b-i  can  be  reduced  to  an  equivalent  partial  sun  by  dividing 
ieM  jiN  J  J 


through  by  (Pqj-Rq  -  Kj)  since  this  expression  has  the  same  value  for  all  jtS. 

Case  b.  Assume  that  the  extra  constraint  is  equivalent  to  a  partial 

sum  of  variables  In  a  single  destination  q  (i.e.  the  constraint  Is  of  the  form 

z  x^q=  f  for  TCM  illustrated  in  Table  1).  The  appropriate  linear 
i^T 

combination  for  this  case  can  be  found  in  the  same  manner  as  In  case  a.  After 
the  origin  constraint  multiplier  R]  is  set  equal  to  zero  and  the  destination 
constraint  multipliers  are  determined  using  the  equations  Rj  +  Kj  =  p^j  ,  the 
remaining  Rj  values  can  be  found.  In  every  origin  i  that  does  not  include 
variables  in  the  equivalent  partial  sum  (i.e.  for  ieM-T),  there  exists  a 
unique  value  for  R-j  such  that  the  constraints  R^  +  Kj  =  p-jj  will  be  satisfied 
simultaneously  by  a  unique  value  for  all  j,;N.  For  those  origins  ieT  that  have 
a  variable  in  common  with  the  equivalent  extra  constraint,  the  equation  R-j  *  Kj=  Pij 
will  be  satisfied  by  a  unique  value  for  all  j  j*  q.  Setting  Ri,  IeT  equal 
to  this  unique  value  will  provide  the  scalar  multipliers  for  a  linear  com¬ 
bination  of  the  node  constraints  which  when  subtracted  from  the  original  extra 
constraint  yield  an  equivalent  partial  sum  of  variables  in  the  single  destination 
constraint  q  associated  with  the  origins  in  T.  The  equivalent  constraint 

z  (Pin-  Ri  -  ^q)  Mcr  d  -  £  Ria-j-  z  Kj bj  can  be  reduced  to  an  equivalent 
ieT  ieM  jf;N 

partial  sum  by  dividing  through  by  (p-jq-R-j  -Kq)  since  this  expression  has 
the  same  value  for  all  itT. 

Case  c.  Assume  the  extra  constraint  is  equivalent  to  a  partial  sum  of 

variables  in  origin  constraint  1  (i.e.  the  constraint  is  z  xj  f  for  ScN 

jeS 

as  illustrated  in  Table  1).  Starting  again  with  R-^  equal  to  zero,  the  Kj 
values  for  jeN  are  immediately  determined.  Since  the  equivalent  partial  sum 
is  assumed  to  be  in  origin  1,  multipliers  must  be  found  to  satisfy  R-j  +Kj  =  pjj 
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1 


for  all  1  ana  j  such  that  ifM  in  order  for  these  variables  not  to  appear  in  the 
equivalent  constraint.  However  given  the  current  values  for  Kj  unique  values 
for  Rj,  1  t  1  cannot  be  found,  for  if  they  could  be  then  Rj  +  Kj  =  p^j  would 
hold  for  all  i  and  all  j  contradicting  the  assumption.  Thus  there  must  be  two 
possible  values  for  each  Rj,  i  M.  Set  these  R^  equal  to  the  value  that 
satisfies  R^  +  Kj  =  p^j  for  Thus,  R^  +  K j ^ P i j  f°r  i  t  1  and  for  jeS. 

Note,  however,  new  values  for  K.  can  be  found  that  satisfy  R.  +  K.  =  p..  for 

J  ■  J  '  J 

i  f  1  and  jcS.  Thus  by  changing  the  appropriate  K.  values  after  having 

J 

determined  values  for  all  R^,  one  can  find  the  linear  combination  of  the 
node  constraints  which  when  subtracted  from  the  extra  constraint  leaves  the 
equivalent  restriction  l  (p^-K.)  xH=  d  -  i  R.A<  -  i  K^b^ . Since ( p-.  .-K^ )  is 

jeS  1 j  J  'J  icM  1  j£N  J  J  J 

the  same  for  all  j*S  then  dividing  through  by  this  quantity  will  yield  the 
equivalent  partial  sum. 

Based  on  the  reasoning  just  presented  we  can  describe  a  general  procedure 
to  effect  the  transformation  as  required  in  the  three  cases.  It  is  convenient 
to  use  a  transportation  tableau  with  the  coefficients  of  the  extra  constraint 
placed  In  the  cell  corresponding  to  the  appropriate  variable.  This  format  is 
Illustrated  in  Table  2.  Values  for  the  R^  and  K.  multipliers  are  shown  around 
the  rim  of  the  tableau  and  can  be  determined  in  a  manner  similar  to  that  used 
to  find  values  for  dual  evaluators  for  a  basic  solution  to  a  transportation 
problem. 


Table  2.  Tableau  Representation 
of  an  Extra  Constraint 


Our  procedure  for  finding  the  appropriate  linear  oombination  of  the 
constraints  is: 

Step  1.  Set  Rj.  equal  to  zero. 

Step  2.  For  all  jrN  set  Kj  equal  to  p^j. 

Step  3.  Try  to  determine  a  unique  value  for  R2  using  the  equations 
r2  +  Kj  =  ?2j  f°r  3eN’ 

a)  If  a  unique  value  for  R2  can  be  found  set  R2  equal  to  this 
value  and  proceed  to  step  4. 

b)  If  the  equations  are  not  satisfied  by  a  unique  value  of  R2  but 
all  equations  except  q  can  be  satisfied  by  a  single  value  of  R2, 

set  R2  so  that  the  ^2  +  Kq  t  P2q  an^  mark  cell  (2, q) with  a  "star". How¬ 
ever,  if  this  is  not  the  case,  and  R2  must  assume  two  distinct  values  to 
satisfy  the  equations  for  all  j,  set  R2  equal  to  either  of  two 
values  arbitrarily  and  "star"  the  cells  for  which  R2+Kj^  P2j*  Proceed 
to  step  4. 

c)  If  more  than  two  values  of  R2  are  required  to  satisfy  all  the 
equations,  stop.  The  constraint  is  not  equivalent  to  a  partial  sun 
of  variables  in  a  single  node  constraint. 

Step  4.  Continue  determining  values  for  the  retaining  R^  as  in  step  2  except 
when  R^  must  assume  two  distinct  values  in  order  for  the  equations 
Ri  +  Kj  =  Pij  to  be  satisfied  for  all  j.  In  this  case  set  R^  so 
that  the  starred  cells  in  this  row  lie  in  the  sane  columns  as  those 
starred  in  earlier  rows.  If  this  cannot  be  done,  stop.  The  eonttlmini 
is  not  equivalent  to  a  partial  sum  of  variables  in  a  single  node 
constraint.  Also  for  any  column  r  with  starred  cells  check  to  see 
that  p^-R^  =  Pkr-Rk  for  iA  where  i  and  k  are  rows  containing  the 
starred  cells.  If  this  is  not  the  case  for  all  ifk,  stop.  Again 
the  constraint  is  not  equivalent  to  a  partial  sum  of  variables  in  a 
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single  node  constraint . 

Step  5.  After  all  Rj  Jutvo  lxxai  determined,  four  cases  are  possible, 
i)  the  starr<d  evils  occur  only  in  a  single  row. 
li)  the  starred  cc.'ls  occur  in  all  cells  in  a  subnet  of  the 
coIutwf  except  tor  the  cells  in  row  1. 
iii)  the  starred  cells  occur  only  in  a  single  column, 
iv)  the  starred  cells  occur  in  some  cells  in  a  subset  of  the 
colunns  but  not  in  row  1  and  not  in  seme  other  rcw. 

In  case  i  the  starred  cells  indicate  the  variables  included  in  the  equi¬ 
valent  partial  sun  of  variables.  For  these  starred  cells  the  coefficient  on 
these  variables  in  an  erjuivalent  constraint  is  p^j-  -It.  After  forming  the 
linear  combination  from  the  extra  constraint  only  those  terms  will  remain.  Since 
pij  ”Ri  “Kj  130  samG  ^or  a11  starred  cells,  the  equivalent  partial  sum 
can  be  obtained  by  dividing  through  by  the  coefficient  on  the  variables. 

In  case  ii  the  Kj  values  for  the  columns  containing  starred  cells  can  be 
changed  so  that  the  aquation  R^  +  K.  =  p-  holds  for  the  starred  cells.  The  effect 
is  to  "erase"  the  stars  from  these  cells  and  place  them  in  the  cells  in  row  1 
in  these  colunns.  This  is  now  the  same  as  case  i. 

In  case  iii  the  starred  cells  indicate  the  variables  in  a  single  destination 
constraint  that  oanpriso  an  equivalent  constraint.  Since  pxj  -R^  -  Kj  will  be 
the  same  for  all  starred  cells  the  equivalent  partial  sun  can  be  obtained  by 
dividing  through  by  the  coefficient  of  the  variables. 

In  case  iv  the  constraint  is  not  equivalent  to  a  partial  sun  of  variables 
in  a  single  origin  or  destination  constraint. 

For  cases  i,  ii,  and  iii  we  have  found  the  coefficients  of  the  variables 

in  an  equivalent  constraint.  The  new  right-hand  side  value  can  be  found  by 

subtracting  j;  +  z  Kjbj  from  the  original  right-hand  side. 

ieM  jeN 
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We  can  summarize  our  results  to  this  [joint  in  tlie  following  tlieoran. 

Theorem:  If  an  extra  constraint  is  equivalent  by  a  linear  transformation 
to  a  partial  sum  of  variables  in  a  single  node  constraint  then  the  stated  pro¬ 
cedure  finds  the  equivalent  partial  sum. 

Proof;  The  hypothesis  allows  us  to  assert  that  the  equivalence  can  be 
determined  by  a  linear  combination  of  the  node  constraints  of  the  transpor¬ 
tation  problem  The  remainder  of  the  proof  is  contained  in  cases  a,b,  and  c  above. 

A  four  origin,  five  destination  constrained  transportation  problem  is  shown 
in  Table  3  with  die  coefficients  of  the  extra  constraint  indicated  in  the  cells 
corresponding  to  the  appropriate  variables  The  associated  supply  and  demand 
values  are  shown  along  the  rim  of  the  table.  Assume  that  the  extra  constraint 
is  an  "equality"  constraint  wi  th  a  right-hand  side  value  of  84. 

Supply 
10 

15 

10 

10 

Demand  8  7  9  6  15 

Table  3 

By  applying  the  procedure  described  above,  tie  following  R^  and  Kj  mul¬ 
tipliers  are  obtained. 

In  step  1  is  set  equal  to  zero. 

In  step  2  the  Kj  multipliers  are  set  equal  to  the  following  values 

Kx  ■  0,  K2  =  1,  K3  =  0,  K4  =  4,  K5  =  -5 

In  step  3  a  unique  value  for  R2  cannot  be  found  to  satisfy  R2  +  Kj  =  p2j 
for  all  j.  The  value  6  satisfies  cfat  equation  for  p2}_  and  P22  and  the  value 
4  satisfies  ttet  equation  for  p23»  P24  and  P25.  Arbitrarily  set  R2  equal  to 
6  and  star  the  cells  p23.  p24,  and  p2,.. 


0 

1 

0 

4 

-5 

6 

7 

4 

8 

-1 

5 

6 

5 

9 

0 

-1 

0 

cLj 

3 

-6 
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As  required  by  step  4  R3  is  set  equal  to  5  because  this  value  satisfies 
the  equation  R3  +  K j  *  P3j  for  all  j.  Similarly  a  value  of  -1  satisfies 
R4  f  Kj  “  P4j  for  all  j.  The  results  of  the  procedure  to  this  point  are  shown 
in  Table  4. 


Ki  -  0  i<2  *  1  K3  =  0  K4  =  4  K5  =  -5  Supply 

Rl  -  0 
R2  *  6 
R3  a  5 
R4  —1 

Demand  8  7  9  6  15 

Table  4 

To  complete  the  transformation  next  compute  the  values  of  (p^  -R^  -Kj ) 
for  the  starred  cells,  check  that  they  are  equal,  and  obtain  the  coefficients 
for  the  variables  in  an  equivalent  constraint.  In  particular  the  left-hand 
side  of  the  equivalent  constraint  is  -2x23  +  -2x24  +  -2x25.  The  new  right- 

4  5 

hand  side  value  is  found  by  the  formula  (+84)  -  z  Ra^-  Z  ICb;  =  84  -  (130)- 

i=l  X  j-1  J  J 

(-44)  =  -2.  Dividing  through  by  -2,  we  obtain  the  equivalent  partial  sum 
x23+x24+x25=  'nius»  ty  using  the  procedure  described  above  we  have  found 
a  partial  sun  of  variables  in  a  single  row  equivalent  to  the  original  extra 
constraint.  The  original  transportation  problan  can  be  enlarged  by  one 
scurce  and  one  destination  in  the  manner  suggested  by  Wagner  [9]  ,and  an  optimal 
solution  to  the  original  problan  can  be  found  using  the  transportation  algorithm* 
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0 

4 

-5  | 

6 

7 

4* 

8* 

-1* 

5 

6 

5 

9 

0 

-1 

0 

-1 

3 

-6 

10 

15 

10 

10 


3.0  BCTENSICNS 

In  the  development  of  the  procedure  for  transforming  extra  constraints  and 
in  the  example  the  original  extra  constraint  was  assumed  to  be  an  "equality" 

type.  It  should  be  clear  that  the  same  transformation  can  be  made  for  both 
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"less  than  or  equal"  and  "greater  than  or  equal"  constraints.  If  there  are 
several  extra  constraints  then  the  procedure  can  be  applied  to  each  one  separ¬ 
ately  to  obtain  an  equivalent  partial  sun  for  each  extra  constraint.  Wagner 
[9]  has  shown  that  if  these  partial  sums  involve  disjoint  sets  of  variables 
and  if  the  sets  are  nested  in  the  same  node  constraints  then  the  problem  can 
be  transformed  into  an  enlarged  transportation  problem. 

4.0  APPLICATIONS 


Many  models  have  the  structure  of  a  transportation  problem  with  additional 
restrictions.  The  extra  constraints  may  represent  secondary  objectives  or 
restrictions  that  are  not  reflected  in  trie  standard  node  constraints.  To 
illustrate  a  typical  example  of  this  class  of  problems,  consider  the  trans¬ 
portation  model  where  warehouses  supply  markets  and  the  objective  is  the  stan¬ 
dard  one  of  finding  a  shipping  pattern  that  will  minimize  the  total  shipping 
oost.  Suppose  additionally  that  the  products  shipped  are  of  a  fragile  nature 
and  if  they  are  sent  via  particular  routes  each  item  mast  be  specially  pack¬ 
aged  to  prevent  losses  in  shipping.  Table  6  shows  the  packaging  time  in 
minutes  per  unit  required  to  prepare  a  unit  for  sliipment  on  the  various  routes. 
Suppose  that  we  wish  to  limit  the  average  packaging  time  per  unit  to  at  most 
2-1/2  minutes  per  unit  (i.e. ,  for  the  60  units  that  must  be  shipped,  the  total 
packaging  time  must  not  exceed  150  minutes).  Using  our  procedure  a  linear 
combination  of  the  node  constraints  of  the  transportation  problems  can  be 
found  which  when  subtracted  from  the  original  extra  constraint  yields  an 
equivalent  partial  sun  of  variables  x22+  *24  i  5*  The  node  constraint 
multipliers  for  the  linear  oanbination  are  indicated  in  table  6.  Thus  the 
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Kl  =  0  ^=1  =  0  =  2  Kg  =  l  Supply 


o 

II 
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0 
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1 

10 

r2  -  i 

1 

3* 

1 

4* 

2 

20 

R3  =  3 

3 

4 

3 

5 

4 

20 

n 

to 

2 

3 

2 

4 

3 

10 

Demand 

15 

10 

5 

5 

25 

60 

Total 


Table  6 


2-1/2  minute  average  packaging  time  restriction  can  only  be  satisfied  if  the 
total  number  of  units  shipped  along  routes  (2,2)  and  (2,4)  is  less  than  or 
equal  to  5.  The  transportation  problem  can  be  transformed  to  include  this 
restriction  directly, and  thus  the  original  constrained  transportation  problem 
can  be  solved  as  a  transportation  problem  witn  one  additional  origin  and  one 
additional  destination. 
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